3.21.93 \(\int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx\) [2093]

3.21.93.1 Optimal result
3.21.93.2 Mathematica [A] (verified)
3.21.93.3 Rubi [A] (verified)
3.21.93.4 Maple [C] (warning: unable to verify)
3.21.93.5 Fricas [A] (verification not implemented)
3.21.93.6 Sympy [F(-1)]
3.21.93.7 Maxima [A] (verification not implemented)
3.21.93.8 Giac [B] (verification not implemented)
3.21.93.9 Mupad [F(-1)]

3.21.93.1 Optimal result

Integrand size = 35, antiderivative size = 152 \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 (b d-a e)^2 (d+e x)^{7/2} \sqrt {a^2+2 a b x+b^2 x^2}}{7 e^3 (a+b x)}-\frac {4 b (b d-a e) (d+e x)^{9/2} \sqrt {a^2+2 a b x+b^2 x^2}}{9 e^3 (a+b x)}+\frac {2 b^2 (d+e x)^{11/2} \sqrt {a^2+2 a b x+b^2 x^2}}{11 e^3 (a+b x)} \]

output
2/7*(-a*e+b*d)^2*(e*x+d)^(7/2)*((b*x+a)^2)^(1/2)/e^3/(b*x+a)-4/9*b*(-a*e+b 
*d)*(e*x+d)^(9/2)*((b*x+a)^2)^(1/2)/e^3/(b*x+a)+2/11*b^2*(e*x+d)^(11/2)*(( 
b*x+a)^2)^(1/2)/e^3/(b*x+a)
 
3.21.93.2 Mathematica [A] (verified)

Time = 0.03 (sec) , antiderivative size = 79, normalized size of antiderivative = 0.52 \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 \sqrt {(a+b x)^2} (d+e x)^{7/2} \left (99 a^2 e^2+22 a b e (-2 d+7 e x)+b^2 \left (8 d^2-28 d e x+63 e^2 x^2\right )\right )}{693 e^3 (a+b x)} \]

input
Integrate[(a + b*x)*(d + e*x)^(5/2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2],x]
 
output
(2*Sqrt[(a + b*x)^2]*(d + e*x)^(7/2)*(99*a^2*e^2 + 22*a*b*e*(-2*d + 7*e*x) 
 + b^2*(8*d^2 - 28*d*e*x + 63*e^2*x^2)))/(693*e^3*(a + b*x))
 
3.21.93.3 Rubi [A] (verified)

Time = 0.25 (sec) , antiderivative size = 99, normalized size of antiderivative = 0.65, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.114, Rules used = {1187, 27, 53, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (a+b x) \sqrt {a^2+2 a b x+b^2 x^2} (d+e x)^{5/2} \, dx\)

\(\Big \downarrow \) 1187

\(\displaystyle \frac {\sqrt {a^2+2 a b x+b^2 x^2} \int b (a+b x)^2 (d+e x)^{5/2}dx}{b (a+b x)}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\sqrt {a^2+2 a b x+b^2 x^2} \int (a+b x)^2 (d+e x)^{5/2}dx}{a+b x}\)

\(\Big \downarrow \) 53

\(\displaystyle \frac {\sqrt {a^2+2 a b x+b^2 x^2} \int \left (\frac {b^2 (d+e x)^{9/2}}{e^2}-\frac {2 b (b d-a e) (d+e x)^{7/2}}{e^2}+\frac {(a e-b d)^2 (d+e x)^{5/2}}{e^2}\right )dx}{a+b x}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {\sqrt {a^2+2 a b x+b^2 x^2} \left (-\frac {4 b (d+e x)^{9/2} (b d-a e)}{9 e^3}+\frac {2 (d+e x)^{7/2} (b d-a e)^2}{7 e^3}+\frac {2 b^2 (d+e x)^{11/2}}{11 e^3}\right )}{a+b x}\)

input
Int[(a + b*x)*(d + e*x)^(5/2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2],x]
 
output
(Sqrt[a^2 + 2*a*b*x + b^2*x^2]*((2*(b*d - a*e)^2*(d + e*x)^(7/2))/(7*e^3) 
- (4*b*(b*d - a*e)*(d + e*x)^(9/2))/(9*e^3) + (2*b^2*(d + e*x)^(11/2))/(11 
*e^3)))/(a + b*x)
 

3.21.93.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 53
Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int 
[ExpandIntegrand[(a + b*x)^m*(c + d*x)^n, x], x] /; FreeQ[{a, b, c, d, n}, 
x] && IGtQ[m, 0] && ( !IntegerQ[n] || (EqQ[c, 0] && LeQ[7*m + 4*n + 4, 0]) 
|| LtQ[9*m + 5*(n + 1), 0] || GtQ[m + n + 2, 0])
 

rule 1187
Int[((d_.) + (e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))^(n_.)*((a_) + (b_.)*(x_ 
) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(a + b*x + c*x^2)^FracPart[p]/(c^ 
IntPart[p]*(b/2 + c*x)^(2*FracPart[p]))   Int[(d + e*x)^m*(f + g*x)^n*(b/2 
+ c*x)^(2*p), x], x] /; FreeQ[{a, b, c, d, e, f, g, m, n, p}, x] && EqQ[b^2 
 - 4*a*c, 0] &&  !IntegerQ[p]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.21.93.4 Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 2.

Time = 0.25 (sec) , antiderivative size = 69, normalized size of antiderivative = 0.45

method result size
default \(\frac {2 \,\operatorname {csgn}\left (b x +a \right ) \left (e x +d \right )^{\frac {7}{2}} \left (63 b^{2} e^{2} x^{2}+154 a b \,e^{2} x -28 b^{2} d e x +99 e^{2} a^{2}-44 a b d e +8 b^{2} d^{2}\right )}{693 e^{3}}\) \(69\)
gosper \(\frac {2 \left (e x +d \right )^{\frac {7}{2}} \left (63 b^{2} e^{2} x^{2}+154 a b \,e^{2} x -28 b^{2} d e x +99 e^{2} a^{2}-44 a b d e +8 b^{2} d^{2}\right ) \sqrt {\left (b x +a \right )^{2}}}{693 e^{3} \left (b x +a \right )}\) \(79\)
risch \(\frac {2 \sqrt {\left (b x +a \right )^{2}}\, \left (63 b^{2} e^{5} x^{5}+154 a b \,e^{5} x^{4}+161 b^{2} d \,e^{4} x^{4}+99 a^{2} e^{5} x^{3}+418 a b d \,e^{4} x^{3}+113 b^{2} d^{2} e^{3} x^{3}+297 a^{2} d \,e^{4} x^{2}+330 a b \,d^{2} e^{3} x^{2}+3 b^{2} d^{3} e^{2} x^{2}+297 a^{2} d^{2} e^{3} x +22 a b \,d^{3} e^{2} x -4 b^{2} d^{4} e x +99 a^{2} d^{3} e^{2}-44 a b \,d^{4} e +8 b^{2} d^{5}\right ) \sqrt {e x +d}}{693 \left (b x +a \right ) e^{3}}\) \(198\)

input
int((b*x+a)*(e*x+d)^(5/2)*((b*x+a)^2)^(1/2),x,method=_RETURNVERBOSE)
 
output
2/693*csgn(b*x+a)*(e*x+d)^(7/2)*(63*b^2*e^2*x^2+154*a*b*e^2*x-28*b^2*d*e*x 
+99*a^2*e^2-44*a*b*d*e+8*b^2*d^2)/e^3
 
3.21.93.5 Fricas [A] (verification not implemented)

Time = 0.34 (sec) , antiderivative size = 174, normalized size of antiderivative = 1.14 \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 \, {\left (63 \, b^{2} e^{5} x^{5} + 8 \, b^{2} d^{5} - 44 \, a b d^{4} e + 99 \, a^{2} d^{3} e^{2} + 7 \, {\left (23 \, b^{2} d e^{4} + 22 \, a b e^{5}\right )} x^{4} + {\left (113 \, b^{2} d^{2} e^{3} + 418 \, a b d e^{4} + 99 \, a^{2} e^{5}\right )} x^{3} + 3 \, {\left (b^{2} d^{3} e^{2} + 110 \, a b d^{2} e^{3} + 99 \, a^{2} d e^{4}\right )} x^{2} - {\left (4 \, b^{2} d^{4} e - 22 \, a b d^{3} e^{2} - 297 \, a^{2} d^{2} e^{3}\right )} x\right )} \sqrt {e x + d}}{693 \, e^{3}} \]

input
integrate((b*x+a)*(e*x+d)^(5/2)*((b*x+a)^2)^(1/2),x, algorithm="fricas")
 
output
2/693*(63*b^2*e^5*x^5 + 8*b^2*d^5 - 44*a*b*d^4*e + 99*a^2*d^3*e^2 + 7*(23* 
b^2*d*e^4 + 22*a*b*e^5)*x^4 + (113*b^2*d^2*e^3 + 418*a*b*d*e^4 + 99*a^2*e^ 
5)*x^3 + 3*(b^2*d^3*e^2 + 110*a*b*d^2*e^3 + 99*a^2*d*e^4)*x^2 - (4*b^2*d^4 
*e - 22*a*b*d^3*e^2 - 297*a^2*d^2*e^3)*x)*sqrt(e*x + d)/e^3
 
3.21.93.6 Sympy [F(-1)]

Timed out. \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\text {Timed out} \]

input
integrate((b*x+a)*(e*x+d)**(5/2)*((b*x+a)**2)**(1/2),x)
 
output
Timed out
 
3.21.93.7 Maxima [A] (verification not implemented)

Time = 0.21 (sec) , antiderivative size = 214, normalized size of antiderivative = 1.41 \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 \, {\left (7 \, b e^{4} x^{4} - 2 \, b d^{4} + 9 \, a d^{3} e + {\left (19 \, b d e^{3} + 9 \, a e^{4}\right )} x^{3} + 3 \, {\left (5 \, b d^{2} e^{2} + 9 \, a d e^{3}\right )} x^{2} + {\left (b d^{3} e + 27 \, a d^{2} e^{2}\right )} x\right )} \sqrt {e x + d} a}{63 \, e^{2}} + \frac {2 \, {\left (63 \, b e^{5} x^{5} + 8 \, b d^{5} - 22 \, a d^{4} e + 7 \, {\left (23 \, b d e^{4} + 11 \, a e^{5}\right )} x^{4} + {\left (113 \, b d^{2} e^{3} + 209 \, a d e^{4}\right )} x^{3} + 3 \, {\left (b d^{3} e^{2} + 55 \, a d^{2} e^{3}\right )} x^{2} - {\left (4 \, b d^{4} e - 11 \, a d^{3} e^{2}\right )} x\right )} \sqrt {e x + d} b}{693 \, e^{3}} \]

input
integrate((b*x+a)*(e*x+d)^(5/2)*((b*x+a)^2)^(1/2),x, algorithm="maxima")
 
output
2/63*(7*b*e^4*x^4 - 2*b*d^4 + 9*a*d^3*e + (19*b*d*e^3 + 9*a*e^4)*x^3 + 3*( 
5*b*d^2*e^2 + 9*a*d*e^3)*x^2 + (b*d^3*e + 27*a*d^2*e^2)*x)*sqrt(e*x + d)*a 
/e^2 + 2/693*(63*b*e^5*x^5 + 8*b*d^5 - 22*a*d^4*e + 7*(23*b*d*e^4 + 11*a*e 
^5)*x^4 + (113*b*d^2*e^3 + 209*a*d*e^4)*x^3 + 3*(b*d^3*e^2 + 55*a*d^2*e^3) 
*x^2 - (4*b*d^4*e - 11*a*d^3*e^2)*x)*sqrt(e*x + d)*b/e^3
 
3.21.93.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 630 vs. \(2 (107) = 214\).

Time = 0.28 (sec) , antiderivative size = 630, normalized size of antiderivative = 4.14 \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 \, {\left (3465 \, \sqrt {e x + d} a^{2} d^{3} \mathrm {sgn}\left (b x + a\right ) + 3465 \, {\left ({\left (e x + d\right )}^{\frac {3}{2}} - 3 \, \sqrt {e x + d} d\right )} a^{2} d^{2} \mathrm {sgn}\left (b x + a\right ) + \frac {2310 \, {\left ({\left (e x + d\right )}^{\frac {3}{2}} - 3 \, \sqrt {e x + d} d\right )} a b d^{3} \mathrm {sgn}\left (b x + a\right )}{e} + 693 \, {\left (3 \, {\left (e x + d\right )}^{\frac {5}{2}} - 10 \, {\left (e x + d\right )}^{\frac {3}{2}} d + 15 \, \sqrt {e x + d} d^{2}\right )} a^{2} d \mathrm {sgn}\left (b x + a\right ) + \frac {231 \, {\left (3 \, {\left (e x + d\right )}^{\frac {5}{2}} - 10 \, {\left (e x + d\right )}^{\frac {3}{2}} d + 15 \, \sqrt {e x + d} d^{2}\right )} b^{2} d^{3} \mathrm {sgn}\left (b x + a\right )}{e^{2}} + \frac {1386 \, {\left (3 \, {\left (e x + d\right )}^{\frac {5}{2}} - 10 \, {\left (e x + d\right )}^{\frac {3}{2}} d + 15 \, \sqrt {e x + d} d^{2}\right )} a b d^{2} \mathrm {sgn}\left (b x + a\right )}{e} + 99 \, {\left (5 \, {\left (e x + d\right )}^{\frac {7}{2}} - 21 \, {\left (e x + d\right )}^{\frac {5}{2}} d + 35 \, {\left (e x + d\right )}^{\frac {3}{2}} d^{2} - 35 \, \sqrt {e x + d} d^{3}\right )} a^{2} \mathrm {sgn}\left (b x + a\right ) + \frac {297 \, {\left (5 \, {\left (e x + d\right )}^{\frac {7}{2}} - 21 \, {\left (e x + d\right )}^{\frac {5}{2}} d + 35 \, {\left (e x + d\right )}^{\frac {3}{2}} d^{2} - 35 \, \sqrt {e x + d} d^{3}\right )} b^{2} d^{2} \mathrm {sgn}\left (b x + a\right )}{e^{2}} + \frac {594 \, {\left (5 \, {\left (e x + d\right )}^{\frac {7}{2}} - 21 \, {\left (e x + d\right )}^{\frac {5}{2}} d + 35 \, {\left (e x + d\right )}^{\frac {3}{2}} d^{2} - 35 \, \sqrt {e x + d} d^{3}\right )} a b d \mathrm {sgn}\left (b x + a\right )}{e} + \frac {33 \, {\left (35 \, {\left (e x + d\right )}^{\frac {9}{2}} - 180 \, {\left (e x + d\right )}^{\frac {7}{2}} d + 378 \, {\left (e x + d\right )}^{\frac {5}{2}} d^{2} - 420 \, {\left (e x + d\right )}^{\frac {3}{2}} d^{3} + 315 \, \sqrt {e x + d} d^{4}\right )} b^{2} d \mathrm {sgn}\left (b x + a\right )}{e^{2}} + \frac {22 \, {\left (35 \, {\left (e x + d\right )}^{\frac {9}{2}} - 180 \, {\left (e x + d\right )}^{\frac {7}{2}} d + 378 \, {\left (e x + d\right )}^{\frac {5}{2}} d^{2} - 420 \, {\left (e x + d\right )}^{\frac {3}{2}} d^{3} + 315 \, \sqrt {e x + d} d^{4}\right )} a b \mathrm {sgn}\left (b x + a\right )}{e} + \frac {5 \, {\left (63 \, {\left (e x + d\right )}^{\frac {11}{2}} - 385 \, {\left (e x + d\right )}^{\frac {9}{2}} d + 990 \, {\left (e x + d\right )}^{\frac {7}{2}} d^{2} - 1386 \, {\left (e x + d\right )}^{\frac {5}{2}} d^{3} + 1155 \, {\left (e x + d\right )}^{\frac {3}{2}} d^{4} - 693 \, \sqrt {e x + d} d^{5}\right )} b^{2} \mathrm {sgn}\left (b x + a\right )}{e^{2}}\right )}}{3465 \, e} \]

input
integrate((b*x+a)*(e*x+d)^(5/2)*((b*x+a)^2)^(1/2),x, algorithm="giac")
 
output
2/3465*(3465*sqrt(e*x + d)*a^2*d^3*sgn(b*x + a) + 3465*((e*x + d)^(3/2) - 
3*sqrt(e*x + d)*d)*a^2*d^2*sgn(b*x + a) + 2310*((e*x + d)^(3/2) - 3*sqrt(e 
*x + d)*d)*a*b*d^3*sgn(b*x + a)/e + 693*(3*(e*x + d)^(5/2) - 10*(e*x + d)^ 
(3/2)*d + 15*sqrt(e*x + d)*d^2)*a^2*d*sgn(b*x + a) + 231*(3*(e*x + d)^(5/2 
) - 10*(e*x + d)^(3/2)*d + 15*sqrt(e*x + d)*d^2)*b^2*d^3*sgn(b*x + a)/e^2 
+ 1386*(3*(e*x + d)^(5/2) - 10*(e*x + d)^(3/2)*d + 15*sqrt(e*x + d)*d^2)*a 
*b*d^2*sgn(b*x + a)/e + 99*(5*(e*x + d)^(7/2) - 21*(e*x + d)^(5/2)*d + 35* 
(e*x + d)^(3/2)*d^2 - 35*sqrt(e*x + d)*d^3)*a^2*sgn(b*x + a) + 297*(5*(e*x 
 + d)^(7/2) - 21*(e*x + d)^(5/2)*d + 35*(e*x + d)^(3/2)*d^2 - 35*sqrt(e*x 
+ d)*d^3)*b^2*d^2*sgn(b*x + a)/e^2 + 594*(5*(e*x + d)^(7/2) - 21*(e*x + d) 
^(5/2)*d + 35*(e*x + d)^(3/2)*d^2 - 35*sqrt(e*x + d)*d^3)*a*b*d*sgn(b*x + 
a)/e + 33*(35*(e*x + d)^(9/2) - 180*(e*x + d)^(7/2)*d + 378*(e*x + d)^(5/2 
)*d^2 - 420*(e*x + d)^(3/2)*d^3 + 315*sqrt(e*x + d)*d^4)*b^2*d*sgn(b*x + a 
)/e^2 + 22*(35*(e*x + d)^(9/2) - 180*(e*x + d)^(7/2)*d + 378*(e*x + d)^(5/ 
2)*d^2 - 420*(e*x + d)^(3/2)*d^3 + 315*sqrt(e*x + d)*d^4)*a*b*sgn(b*x + a) 
/e + 5*(63*(e*x + d)^(11/2) - 385*(e*x + d)^(9/2)*d + 990*(e*x + d)^(7/2)* 
d^2 - 1386*(e*x + d)^(5/2)*d^3 + 1155*(e*x + d)^(3/2)*d^4 - 693*sqrt(e*x + 
 d)*d^5)*b^2*sgn(b*x + a)/e^2)/e
 
3.21.93.9 Mupad [F(-1)]

Timed out. \[ \int (a+b x) (d+e x)^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\int \sqrt {{\left (a+b\,x\right )}^2}\,\left (a+b\,x\right )\,{\left (d+e\,x\right )}^{5/2} \,d x \]

input
int(((a + b*x)^2)^(1/2)*(a + b*x)*(d + e*x)^(5/2),x)
 
output
int(((a + b*x)^2)^(1/2)*(a + b*x)*(d + e*x)^(5/2), x)